Photoionization of helium-like ions in asymptotic nonrelativistic region 
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Abstract 

The cross section for single K-shell ionization by a high-energy photon is evaluated in the next-to-leading 
order of the nonrelativistic perturbation theory with respect to the electron-electron interaction. The screen- 
ing corrections are of particular importance for light helium-like ions. Even in the case of neutral He atom, 
our analytical predictions turn out to be in good agreement with the numerical calculations performed with 
the use of the sophisticated wave functions. The universal high-energy behavior is studied for the ratio of 
double-to-single photoionization cross sections. We also discuss the fast convergence of the perturbation 
theory over the reversed nuclear charge number 1/Z. 

PACS numbers: 32.80.Fb, 32.80.-t, 31.25.Eb 
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The single and double photoeffects on the helium isoelectronic sequence represent the simplest 
fundamental processes, which are being intensively investigated during last decades 0, ls,0 0,0,0] • 
The accurate treatment of electron correlations is still one of the main theoretical problems. Due 
to recent developments of novel synchrotron radiation sources, the study of the ionization of inner- 
shell electrons by high-energy photons is of particular interest. 

For the theoretical description of the ionization processes on light atomic systems, it is usual to 
employ sophisticated methods with highly correlated wave functions. This allows one to take into 
account electron correlation effects beyond the independent-particle approximation. However, all 
the methods suffer from the gauge dependence. The latter can serve as a level of accuracy for the 
theoretical predictions. In addition, the final results for cross sections of the ionization processes 
are presented in a numerical form, which is not always easy to analyze. 

In the case of heavy multicharged ions, on the contrary, the usual starting point is the approx- 
imation of non- interacting electrons, which are described by the Coulomb wave functions for the 
discrete and continuous spectra. The electron-electron interaction is treated within the framework 
of perturbation theory, which is also referred to as the expansion with respect to the parameter 
\jZ. The latter represents the ratio of the strength of the electron-electron interaction to the 
electron-nucleus one. To leading orders, perturbation theory allows one to derive analytical re- 
sults. Accounting for higher-order correlation corrections improves the accuracy of the analytical 
predictions in the domain of lower values of the nuclear charge number Z. The results obtained 
within the framework of perturbation theory for the binding energies and for the cross sections are 
gauge independent. 

In this Letter, we evaluate the next-to-leading-order correlation correction to the cross section 
for single K-shell ionization at asymptotic photon energies u characterized by / <C to <C m, where 
I = if /(2m) is the Coulomb potential for single ionization with rj = maZ being the average 
momentum of a K-shell electron, m is the electron mass, and a is the fine-structure constant 
(h = 1, c = 1). The ejected electrons are considered as being nonrelativistic. Accordingly, the 
Coulomb parameter is supposed to be sufficiently small, that is, aZ -C 1. 

Neglecting terms of order (aZ) 2 , the operator describing the electron-photon interaction reads 
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Here p is the momentum operator of an electron, which, in the coordinate representation, is cast into 
the gradient form p = — iV. An incoming photon is characterized by the momentum k, the energy 
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to = \k\ = k, and the polarization vector e. We employ the Coulomb gauge, in which (e • k) = 
and (e* • e) = 1. In general, the nonrelativistic interaction between an electron and a photon 
includes also spin-dependent terms. However, in the case of single and double photoeffects the 
corresponding contributions to the cross sections are strongly suppressed 0, 0, Q] and, therefore, 
can be neglected. 

In the nonrelativistic approximation, spatial and spin parts of two-electron wave functions 
factorize. Moreover, the operator (J2) does not involve the spin matrices. As a result, the spin 
functions can be omitted throughout this consideration, while the symmetry of the coordinate 
wave functions \&(ri,r2) is preserved in the ionization process. The total amplitude of the single 
photoeffect on helium-like ion is given by 

A = 2(V f \V Y \V i ), (2) 

where the factor 2 accounts for the one-particle character of the operator That is, it is sufficient 
to consider the interaction of an incoming photon with a single atomic electron only. 

In first-order perturbation theory with respect to the electron-electron interaction, the wave 
functions of the initial and final states are represented as follows ~ tyf'J + ^^j- Ac- 
cordingly, the amplitude of the process is just A — A^ + .A'- 1 ). Neglecting the electron-electron 
interaction, we shall employ the single-particle approximation in the external Coulomb field of the 
nucleus (Furry picture): 

*j 0) (ri,r 2 ) = Vi.(ri)^,(r 2 ), (3) 
tfJVi.ra) = -J= [V> P (n)Vi s (r 2 ) + M r 2)Mri)} . (4) 



Here p is the momentum of escaping electron at infinity. The explicit expressions for the first-order 




corrections Vf- 1 ) to the wave functions can be found in works [1CL 111!. 
In zeroth approximation, the amplitude (0) looks as follows 

A® = 2(tff Ifytff) = V2{iP p \V^ ls ) . (5) 

The matrix element © can be represented by the Feynman graph depicted in Fig. ^a). Apart 
from the common factor, the expression (JSJ) coincides with the amplitude for single photoeffect on 
hydrogen-like ion in the ground state. 

In the asymptotic region, the energy and momentum conservation laws for the single K-shell 
photoionization keep the same form for both H- and He-like ions: 

E p = E u + uj~uj, (6) 
p = q + k. (7) 



Here E\ s = — I is the Coulomb energy of the K-shell electron and — q is the recoil momentum 
transferred to the nucleus. 

Let us evaluate the amplitude © using the momentum representation. Within the Born ap- 
proximation, the wave function of the ejected high-energy electron is described by a plane wave, 
that is, 

W P )^</|P> = (27T)^(/-P)- (8) 

Here the standard normalization on 5 function in the momenta is employed. Then the matrix 
element © yields 

A^ = V2N y (q\ihs)(e-p). (9) 

If q S> rj holds, the Coulomb wave function of the K-shell electron reads 

(q\^ u )^N ls ^, (10) 

where Nf s = r] 5 /ir. Since in the nonrelativistic domain the momentum k of a photon is negligibly 
small with respect to the electron momentum p, relation Q can be written as 

(11) 

The latter is equivalent to the use of the dipole approximation. Then one can set q 2 = p 2 = 2mu> 
into Eq. (|l(Jj) . The amplitude @ can be further simplified. The corresponding expression for the 
total cross section is well known [3] 

2 8 vtq m 5/2 

- , (I«w«m). (12) 
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If the photon energy is not too high, as the wave function of the final state one needs to utilize 
the one-electron Coulomb wave function of the continuous spectrum. Accordingly, within the dipole 
approximation the amplitude (JSJ) leads to the following expression for the total cross section, 

+ 2 2 10 tt 2 exp(-4gcot- 1 Q ( ^ s m , 

^ =aQ ° 3Z 2 4[l-exp(-27rO] ' (t " < "° ' ( } 



which is in fact valid in the whole nonrelativistic domain Here e 7 = to /I denotes the dimen 



sionless energy for the incident photon, ao = l/(ma) is the Bohr radius, and £ = 1/y^e-y — 1. The 
dimensionless parameter £ _1 = p/rj has the meaning of the momentum p of the ejected electron, 
which is calibrated in units of the characteristic momentum rj. Formula (|12j) provides just the lead- 
ing term in the expansion of Eq. ()13|) with respect to the parameter £ <C 1. Since the expression 



f|13|) involves the combination 7r£, the convergence of the £ expansion is slow. Note also that both 
cross sections (|12jl and are twice as large as that for the single photoeffect on a hydrogen-like 
ion in the ground state. This is the consequence of the approximation of non-interacting electrons 
employed in the derivation. Neither nor account for electron correlation effects. 

Now we shall consider in more details the evaluation of the next-to-leading-order correction .A^ 1 ) 
to the amplitude of single K-shell photoeffect on helium-like ion. In the high-energy nonrelativistic 
limit, the dominant contribution to the amplitude of the process arises only from the Feynman 
diagram depicted in Fig. ^b), providing the Coulomb gauge is employed. This graph accounts 
for the interaction between the electrons in the initial state. All other diagrams, namely, the one, 
which accounts for the electron-electron interaction in the final state, together with both exchange 
diagrams, turn out to be suppressed by the factor of about I/uj and, therefore, can be neglected. 
Accordingly, we can write 

A w =V2{^ ls \V 1 G R (E ls )V 12 \^ ls ) . (14) 

Here the operator V12 describes the Coulomb electron-electron interaction. In the coordinate 
representation, it reads 

V12 = 1 ° 1 • (15) 
\r\ - r 2 \ 

The reduced Green's function G R (Ei s ) corresponding to the energy E\ s of the K-shell electron is 
related to the usual nonrelativistic Coulomb Green's function Gq(E) as follows 

G^) = jta {<*,(*> -M<£1}. o*> 

Within the Born approximation (jHJ, the amplitude ifTljl yields 

= 4W2iV 7 (e-p) A<(|G E (E 1 ,)|/ 1 )i(/i + /|Vi,)x 

where q = p — k. Integrating over the intermediate momenta in Eq. (|17|) . one receives 

AW = A7raV2N y Nf s (e • p) X 

~ ) ^(q\G R (E ls ) (V lv - V i(r]+tl) ) |0), M=2r) , (18) 

(f'\Vix\f) LT 



(/'-/) 2 + A2- 

After taking the derivative with respect to fi, one should set /i = 2rj, where rj = maZ. 



In Eq. (|18[). we shall evaluate first the matrix element with the Coulomb Green's function. Since 



q S> rj, one can employ the integral representation 
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(q\G c (E)Va\0) * 4m V ^ (^Y C (0\V piy+iX \0)dy= (19) 

r ipi f 1 t~ iC -dt , . 

= 2W H [A(l- t )- ipi (l + t )p - < 2 °» 

where pi = \j2mE and £ = 77/pi. In order to isolate the pole contribution together with finite 
terms at E — * E\ s in Eq. ([2()jl. we set i£ = Ei s — 5. The shift 5 is supposed to be small and positive. 
Expanding the intermediate momentum p\ into a series over the parameter 5, one receives 

p a ~ ir/(l + e - e 2 /2) , (21) 
1 — «C ^ e(l-3e/2), (22) 

where e = 5/(21). Performing the integration in Eq. Q2U)) by parts and using the expansions 1)21(1 
and 1)22)1 . we find that 

" + o + TTT ) + ( 23 ) 

(24) 



r? 7 [A(l-t)-ipi(l + t)] 2 (A + r?) 2 Ve 2 A + 77 



+2(A-r?) / 
J 



[A + 7? - (A - ?7)t] 3 ' 

In the integral 1)24)1 we have set 5 = 0, which is equivalent to the substitution ipi — > —77. As a 
result, the matrix element 1)19)1 can be cast into the form 

iq \G C ( Eu - mm = *™£ {-^ Q + § + ^2) + 

1 ln(l-t)eft 
o [277 + (A - 77)tp 



-2(A-r/)/ +0( e )^. (25) 



The analogous matrix element involving the reduced Green's function can be evaluated by 
making use of the definition ()16l) : 

(q\G R (E ls )V iX \0) = VmUq\Gc(Eu-6)Vix\0) + ^(q\^is)(tpu\V lX \0)^ . (26) 

The counter-term in Eq. 1)26)1 is given by 

Here the explicit expression for the matrix element 

47T 

<«fe|0>=JVi.^^ (28) 



and Eq. (fTU|) have been employed. 

Adding Eqs. (|23|) and (|27|). one observes that the pole terms cancel each other. Accordingly, 
we arrive at the following expression 

(q\G R (E ls )V iX \0) = -2 5 vrm^| 2(A ^ )2 + 

' (A_r?) 2(A-H) T ^7^,1 ■ (29) 



(A + ?]) 3 ^ "Jo [27 ? + (A- ?7 )t] 3 
Using Eq. ()29|) allows one to evaluate analytically the matrix element entering Eq. ()18|) . It yields 

Nl (-^) -^(g\G R (E ls ) (V lv - V i{r}+il) ) |0)| M=2 , = ^ai . (30) 

Here the coefFicient ai appears as 

19 S 

ai = -— + -ln2~ -0.6676. (31) 
16 4 

The next-to-leading-order correction A^ to the amplitude of the single photoionization of 
heliumdike ion in the ground state is given by 

=y/2N u N^(e-p)a 1 . (32) 
Lq 

Employing Eqs. Q and l|l()j) yields the total amplitude 

A = A^+A^=A^{l + ^) , (33) 

which accounts for the electron correlations. Within the same approximation, the total cross section 
for the single K-shell photoionization reads 

a + = a+ (I + mZ- 1 ) , (34) 

where cr^ is given by Eq. 1)121) . 

The negative sign of the coefficient (|3Tj) can be understood on the qualitative ground. It is well 
known that the single photoeffect does not proceed on the free electron, while it can occur on the 
bound one Q|. The explanation of this fact follows from the relation (|llj) : the nucleus serves as 
an absorber of the recoil momentum q. In the asymptotic region, the value of q is relatively large, 
since the condition q 3> r] holds. The Coulomb calculation performed within the approximation 
of non-interacting electrons overestimates the cross section, because the electron-nucleus bindings 
are utmost strong in this case. Accounting for the screening effects attenuates these bindings, so 
that, the cross section of the single photoionization reduces in the absolute value. 
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As we already mentioned, if the photon energy is not too high, the cross section crj is more 
preferable rather than . Due to slow convergence of the expansion of Eq. (|13j) with respect to 
the parameter £, it is still legitimate to use of the following formula 

a+ = a+(l + a 1 Z- 1 ) (35) 

instead of Eq. 

As a testing ground for our analytical results, we choose the neutral He atom, which seems to 
be the most thoroughly investigated two-electron system. However, although extensive numeri- 
cal calculations of the photoionization cross sections have been published in the literature, there 
are significant disagreements between predictions based on different sophisticated methods at high 
photon energies. The difficulty of comparison with experimental data arises due to presence of 
additional contributions from the scattering channels. One measures the total attenuation cross 
section, b nt not e„ tne p„o_on one QQ. As tne n,o st _ t_, ca^ 
tions of the photoionization cross sections at high-energy domain, Samson et. al. lj have selected 
results of the work . Bell and Kingston used the Hartree-Fock wave function for the continuum 
state and many-parametrical variational wave function for the ground state. At uj ~ 3.5 keV their 
result for the cross section is equal to 5.74 b, which however exhibits a gauge dependence on the 
level of about 6% |l3j. Our analytical formulas (|34"|) and (|35|) yield 8.24 b and 5.73 b, respectively. 
In this case, the parameter £ = 1/8, but the value 7r£ is not too small. 

In the double photoeffect, one is usually interested in the ratio of double-to-single ionization 
cross sections. At high photon energies, the calculations performed within the framework of the 
leading-order perturbation theory yield 

Ro = -^i- = -7pi > (36) 

where B = 0.090 Q Q and <7q" is given by Eq. (|12|) . Taking into account the higher-order 
screening corrections to the total cross sections leads to the following expression for the universal 
asymptotic ratio 

U (l + a 1 Z- 1 +a 2 Z- 2 + ...) V ; 

The factor Rq is separated out here, since the electron binding energy is supposed to be negligibly 
small compared to the photon energy. To any given order of the perturbation theory with respect 
to the electron-electron interaction the representation (|37|) is the Pade approximant. The next-to- 
leading-order coefficient a\ is given by Eq. (|31|) . while the other coefficients remain to be calculated. 



Nevertheless, employing experimental data for the double-to-single photoionization ratio, one can 
deduce an estimate for the value of the coefficient b\. The latter can be obtained by equating 
the experimental value R exp = 1.72(12)% measured for He atom [3J and the theoretical ratio (|37|) 
truncated with taking into account only the next-to-leading-order correlation corrections. It yields 



fei = -0.981(71) . (38) 

Having fixed the coefficients a\ and bi, we have calculated the double-to-single photoionization 
ratio for helium isoelectronic sequence. In Table [H we present a comparison of our next-to- leading- 
order predictions according to Eq. 1)37(1 with the numerical results obtained by Forrey et al. 4]. The 
account of the screening corrections improves significantly the asymptotic behavior for the leading- 
order ratio Rq i n the case of light two-electron systems. This supports the statement concerning 

n, 

the fast convergence of the 1/Z expansion even in the extreme nonrelativistic domain |l6j]. Indeed, 
the starting approximation of the perturbation theory (assuming non-interacting electrons) would 
be utmost inadequate for the description of light helium-like ions, which are highly correlated. The 
relatively large values for the coefficients ()31j) and (|38[) allow to correct the zeroth approximation. 
Note also that the double-to-single photoionization ratio (|37[l turns out to be less sensitive to the 
higher-order screening corrections rather than the total cross sections, since the coefficients a\ and 
b\ have the same sign. Unfortunately, the significant uncertainty of the coefficient ()38|) distorts 
the true behavior of the ratio R in the case of H ion. The coefficients 02, 61, and 62 should be 
calculated exactly within the framework of the consistent perturbation theory. 

Nevertheless, it is worthwhile to trace out the nontrivial behavior of the series over the parameter 
1/Z taking the binding energy for the ground state in H~ ion as an example. Without the electron- 
electron interaction the Coulomb binding energy is equal to 27.2116 eV. The corrections due to one-, 



ia 



18], and 



two-, and three-photon exchange diagrams are known to yield —17.0079 eV, 4.2921 eV 

—0.1723 eV Q], respectively. Then the total binding energy turns out to be equal to 14.324 eV. 

This should be compared with the exact numerical result of 14.361 eV, which has been obtained 

□ 

within the approximation of an infinitely heavy nucleus [191 ] . On the level of accuracy of about 
10~ 2 eV one already needs to take into account the effect of nuclear recoil. As seen, the terms 
of the 1/Z expansion exhibit sign-changing oscillations and decrease fast in their absolute value. 
Although the formal parameter of the perturbation theory is equal to 1, the actual expansion turns 
out to converge by one order of magnitude due to hidden parameters of the theory. 

Concluding, we have evaluated the single K-shell photoionization cross section with taking into 
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account the next-to-leading-order correlation correction. This allows one to improve the accuracy 
of analytical predictions for light helium-like ions at high-energy domain. We have discussed the 
universal behavior of the double-to-single photoionization ratio as well as the fast convergence of 
the perturbation theory with respect the parameter 1/Z. 
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I 

(a) (b) 

FIG. 1: Fcynman diagrams for the single ionization of K-shell electrons by a single photon. Solid lines 
denote electrons in the Coulomb field of the nucleus, the dashed line denotes the electron-electron Coulomb 
interaction, and the wavy line denotes an incident photon. The line with a heavy dot corresponds to the 
reduced Coulomb Green's function. Diagram (a) does not account for the electron-electron interaction, while 
diagram (b) accounts for it in the initial state. 



TABLE I: The asymptotic ratios of double-to-single photoionization cross sections (in %) are tabulated for 
various nuclear charge numbers Z. The ratio Rq is calculated using the leading-order perturbation theory 
[t! ll"^ . The ratio R is calculated according to the truncated expression (|37l) . taking into account the 
next-to-leading-order corrections with the coefficients (|31() and Q38[l . The numerical calculations by Forrey 
et. al. |4j have been performed with the use of the fully correlated variational wave functions. 
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